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Abstract 

Transverse momentum spectra of charged particles produced in deep inelastic scat- 
tering are measured as a function of the kinematic variables x and Q 2 using the 
HI detector at the ep collider HERA. The data are compared to different parton 
emission models, either with or without ordering of the emissions in transverse mo- 
mentum. The data provide evidence for a relatively large amount of parton radiation 
between the current and the remnant systems. 



HI Collaboration 



C. Adloff 35 , S. Aid 13 , M. Anderson 23 , V. Andreev 26 , B. Andrieu 29 , A. Babaev 25 , 

J. Bahr 36 , J. Ban 18 , Y. Ban 28 , P. Baranov 26 , E. Barrelet 30 , R. Barschke 11 , W. Bartel 11 , 
M. Barth 4 , U. Bassler 30 , H.P. Beck 38 , M. Beck 14 , H.-J. Behrend 11 , A. Belousov 26 , 
Ch. Berger 1 , G. Bernardi 30 , G. Bertrand-Coremans 4 , M. Besangon 9 , R. Beyer 11 , 
P. Biddulph 23 , P. Bispham 23 , J.C. Bizot 28 , V. Blobel 13 , K. Borras 8 , F. Botterweck 27 , 
V. Boudry 29 , A. Braemer 15 , W. Braunschweig 1 , V. Brisson 28 , W. Bruckner 14 , P. Bruel 29 , 

D. Bruncko 18 , C. Brune 16 , R. Buchholz 11 , L. Biingener 13 , J. Biirger 11 , F.W. Biisser 13 , 
A. Buniatian 4 , S. Burke 19 , M.J. Burton 23 , D. Calvet 24 , A.J. Campbell 11 , T. Carli 27 , 
M. Charlet 11 , D. Clarke 5 , A.B. Clegg 19 , B. Clerbaux 4 , S. Cocks 20 , J.C Contreras 8 , 

C. Cormack 20 , J.A. Coughlan 5 , A. Courau 28 , M.-C. Cousinou 24 , G. Cozzika 9 , 

L. Criegee 11 , D.G. Cussans 5 , J. Cvach 31 , S. Dagoret 30 , J.B. Dainton 20 , W.D. Dau 17 , 
K. Daum 42 , M. David 9 , C.L. Davis 19 - 39 , B. Delcourt 28 , A. De Roeck 11 , E.A. De Wolf 4 , 
M. Dirkmann 8 , P. Dixon 19 , P. Di Nezza 33 , W. Dlugosz 7 , C. Dollfus 38 , K.T. Donovan 21 , 
J.D. Dowell 3 , H.B. Dreis 2 , A. Droutskoi 25 , O. Diinger 13 , H. Duhm 12 -* J. Ebert 35 , 
T.R. Ebert 20 , G. Eckerlin 11 , V. Efremenko 25 , S. Egli 38 , R. Eichler 37 , F. Eisele 15 , 

E. Eisenhandler 21 , E. Elsen 11 , M. Erdmann 15 , W. Erdmann 37 , A.B. Fahr 13 , L. Favart 28 , 
A. Fedotov 25 , R. Feist 11 , J. Feltesse 9 , J. Ferencei 18 , F. Ferrarotto 33 , K. Flamm 11 , 

M. Fleischer 8 , M. Flieser 27 , G. Fliigge 2 , A. Fomenko 26 , J. Formanek 32 , J.M. Foster 23 , 
G. Franke 11 , E. Fretwurst 12 , E. Gabathuler 20 , K. Gabathuler 34 , F. Gaede 27 , J. Garvey 3 , 
J. Gayler 11 , M. Gebauer 36 , H. Genzel 1 , R. Gerhards 11 , A. Glazov 36 , L. Goerlich 6 , 
N. Gogitidze 26 , M. Goldberg 30 , D. Goldner 8 , K. Golec-Biernat 6 , B. Gonzalez-Pineiro 30 , 
I. Gorelov 25 , C. Grab 37 , H. Grassier 2 , T. Greenshaw 20 , R.K. Griffiths 21 , 

G. Grindhammer 27 , A. Gruber 27 , C. Gruber 17 , T. Hadig 1 , D. Haidt 11 , L. Hajduk 6 , 
T. Haller 14 , M. Hampel 1 , W.J. Haynes 5 , B. Heinemann 13 , G. Heinzelmann 13 , 
R.C.W. Henderson 19 , H. Henschel 36 , I. Herynek 31 , M.F. Hess 27 , K. Hewitt 3 , 

W. Hildesheim 11 , K.H. Hiller 36 , CD. Hilton 23 , J. Hladky 31 , M. Hoppner 8 , 

D. Hoffmann 11 , T. Holtom 20 , R. Horisberger 34 , V.L. Hudgson 3 , M. Hiitte 8 , 
M. Ibbotson 23 , H. Itterbeck 1 , A. Jacholkowska 28 , C. Jacobsson 22 , M. Jaffre 28 , 

J. Janoth 16 , D.M. Jansen 14 , T. Jansen 11 , L. Jonsson 22 , D.P. Johnson 4 , H. Jung 22 , 
PI.P Kalmus 21 , M. Kander 11 , D. Kant 21 , R. Kaschowitz 2 , U. Kathage 17 , J. Katzy 15 , 

H. H. Kaufmann 36 , O. Kaufmann 15 , M. Kausch 11 , S. Kazarian 11 , I.R. Kenyon 3 , 

S. Kermiche 24 , C. Keuker 1 , C. Kiesling 27 , M. Klein 36 , C. Kleinwort 11 , G. Knies 11 , 
T. Kdhler 1 , J.H. Kohne 27 , H. Kolanoski 36 - 41 , S.D. Kolya 23 , V. Korbel 11 , P. Kostka 36 , 
S.K. Kotelnikov 26 , T. Kramerkamper 8 , MW. Krasny 6 ' 30 , H. Krehbiel 11 , D. Kriicker 27 , 
H. Kiister 22 , M. Kuhlen 27 , T. Kurca 36 , J. Kurzhofer 8 , D. Lacour 30 , B. Laforge 9 , 
M.P.J. Landon 21 , W. Lange 36 , U. Langenegger 37 , A. Lebedev 26 , F. Lehner 11 , 
S. Levonian 29 , G. Lindstrom 12 , M. Lindstroem 22 , F. Linsel 11 , J. Lipinski 13 , B. List 11 , 
G. Lobo 28 , P. Loch 11 - 43 , JW. Lomas 23 , G.C. Lopez 12 , V. Lubimov 25 , D. Luke 8 - 11 , 
L. Lytkin 14 , N. Magnussen 35 , E. Malinovski 26 , R. Maracek 18 , P. Marage 4 , J. Marks 24 , 
R. Marshall 23 , J. Martens 35 , G. Martin 13 , R. Martin 20 , H.-U. Martyn 1 , J. Martyniak 6 , 
T. Mavroidis 21 , S.J. Maxfield 20 , S.J. McMahon 20 , A. Mehta 5 , K. Meier 16 , F. Metlica 14 , 
A. Meyer 13 , A. Meyer 11 , H. Meyer 35 , J. Meyer 11 , P.-O. Meyer 2 , A. Migliori 29 , 
S. Mikocki 6 , D. Milstead 20 , J. Moeck 27 , F. Moreau 29 , J.V. Morris 5 , E. Mroczko 6 , 



f 



D. Miiller 38 , G. Miiller 11 , K. Miiller 11 , P. Murin 18 , V. Nagovizin 25 , R. Nahnhauer 36 , 

B. Naroska 13 , Th. Naumann 36 , I. Negri 24 , P.R. Newman 3 , D. Newton 19 , H.K. Nguyen 30 , 
T.C. Nicholls 3 , F. Niebergall 13 , C. Niebuhr 11 , Ch. Niedzballa 1 , H. Niggli 37 , G. Nowak 6 , 
G.W. Noyes 5 , T. Nunnemann 14 , M. Nyb erg- Wert her 22 , M. Oakden 20 , H. Oberlack 27 , 
J.E. Olsson 11 , D. Ozerov 25 , P. Palmen 2 , E. Panaro 11 , A. Panitch 4 , C. Pascaud 28 , 

G. D. Patel 20 , H. Pawletta 2 , E. Peppel 36 , E. Perez 9 , J.P Phillips 20 , A. Pieuchot 24 , 
D. Pitzl 37 , G. Pope 7 , B. Povh 14 , S. Prell 11 , K. Rabbertz 1 , G. Radel 11 , P. Reimer 31 , 
S. Reinshagen 11 , H. Rick 8 , F. Riepenhausen 2 , S. Riess 13 , E. Rizvi 21 , P. Robmann 38 , 

H. E. Roloff 36,t , R. Roosen 4 , K. Rosenbauer 1 , A. Rostovtsev 25 , F. Rouse 7 , C. Royon 9 , 
K. Riiter 27 , S. Rusakov 26 , K. Rybicki 6 , D.P.C. Sankey 5 , P. Schacht 27 , S. Schiek 13 , 

S. Schleif 16 , P. Schleper 15 , W. von Schlippe 21 , D. Schmidt 35 , G. Schmidt 13 , L. Schoeffel 9 , 
A. Schdning 11 , V. Schroder 11 , E. Schuhmann 27 , B. Schwab 15 , F. Sefkow 38 , R. Sell 11 , 
A. Semenov 25 , V. Shekelyan 11 , I. Sheviakov 26 , L.N. Shtarkov 26 , G. Siegmon 17 , 
U. Siewert 17 , Y. Sirois 29 , I.O. Skillicorn 10 , P. Smirnov 26 , V. Solochenko 25 , Y. Soloviev 26 , 
A. Specka 29 , J. Spiekermann 8 , S. Spielman 29 , H. Spitzer 13 , F. Squinabol 28 , P. Steffen 11 , 
R. Steinberg 2 , H. Steiner 11 ' 40 , J. Steinhart 13 , B. Stella 33 , A. Stellberger 16 , J. Stier 11 , 
J. Stiewe 16 , U. Stofilein 36 , K. Stolze 36 , U. Straumann 15 , W. Struczinski 2 , J.P. Sutton 3 , 
S. Tapprogge 16 , M. Tasevsky 32 , V. Tchernyshov 25 , S. Tchetchelnitski 25 , J. Theissen 2 , 

C. Thiebaux 29 , G. Thompson 21 , N. Tobien 11 , R. Todenhagen 14 , P. Truol 38 , 

G. Tsipolitis 37 , J. Turnau 6 , J. Tutas 15 , E. Tzamariudaki 11 , P. Uelkes 2 , A. Usik 26 , 
S. Valkar 32 , A. Valkarova 32 , C. Vallee 24 , D. Vandenplas 29 , P. Van Esch 4 , 
P. Van Mechelen 4 , Y. Vazdik 26 , P. Verrecchia 9 , G. Villet 9 , K. Wacker 8 , A. Wagener 2 , 
M. Wagener 34 , B. Waugh 23 , G. Weber 13 , M. Weber 16 , D. Wegener 8 , A. Wegner 27 , 
T. Wengler 15 , M. Werner 15 , L.R. West 3 , T. Wilksen 11 , S. Willard 7 , M. Winde 36 , 
G.-G. Winter 11 , C. Wittek 13 , M. Wobisch 2 , E. Wiinsch 11 , J. Zacek 32 , D. Zarbock 12 , 
Z. Zhang 28 , A. Zhokin 25 , P. Zini 30 , F. Zomer 28 , J. Zsembery 9 , K. Zuber 16 , and 
M. zurNedden 38 



1 /. Physikalisches Institut der RWTH, Aachen, Germany 0, 

2 III. Physikalisches Institut der RWTH, Aachen, Germany 0, 

3 School of Physics and Space Research, University of Birmingham, Birmingham, U¥^ 

4 Inter-University Institute for High Energies ULB-VUB, Brussels; Universitaire 
Instelling Antwerpen, Wilrijk; Belgium 

5 Rutherford Appleton Laboratory, Chilton, Didcot, UK** 

6 Institute for Nuclear Physics, Cracow, Poland d 

7 Physics Department and IIRPA, University of California, Davis, California, USA e 

8 Institut fur Physik, Universitdt Dortmund, Dortmund, Germany" 1 

9 CEA, DSM/DAPNIA, CE-Saclay, Gif-sur-Yvette, France 

10 Department of Physics and Astronomy, University of Glasgow, Glasgow, UK* 

11 DESY, Hamburg, Germany 0, 

12 /. Institut fur Experimentalphysik, Universitdt Hamburg, Hamburg, Germany 0, 

13 II. Institut fur Experimentalphysik, Universitdt Hamburg, Hamburg, Germany 0, 

14 Max-Planck-Institut fur Kernphysik, Heidelberg, Germany 

15 Physikalisches Institut, Universitdt Heidelberg, Heidelberg, Germany® 

16 Institut fur Hochenergiephysik, Universitdt Heidelberg, Heidelberg, Germany 0, 



2 



17 Institut fur Reine und Angewandte Kernphysik, Universitdt Kiel, Kiel, Germany " 

18 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak 
Republic?^ 

19 School of Physics and Chemistry, University of Lancaster, Lancaster, UK b 

20 Department of Physics, University of Liverpool, Liverpool, UK° 

21 Queen Mary and Westfield College, London, UK" 

22 Physics Department, University of Lund, Lund, Sweden 9 

23 Physics Department, University of Manchester, Manchester, UL^ 

24 CPPM, Universite d'Aix-Marseille II, IN2P3-CNRS, Marseille, France 

25 Institute for Theoretical and Experimental Physics, Moscow, Russia 

26 Lebedev Physical Institute, Moscow, Russia? 

27 Max-Planck-Institut fur Physik, Miinchen, Germany" 1 

28 LAL, Universite de Paris-Sud, IN2P3-CNRS, Orsay, France 

29 LPNHE, Ecole Poly technique, IN2P3-CNRS, Palaiseau, France 

30 LPNHE, Universites Paris VI and VII, IN2P3-CNRS, Paris, France 

31 Institute of Physics, Czech Academy of Sciences, Praha, Czech Republic?' 11 

32 Nuclear Center, Charles University, Praha, Czech Republic?' 11 

33 INFN Roma 1 and Dipartimento di Fisica, Universita Roma 3, Roma, Italy 

34 Paul Scherrer Institut, Villigen, Switzerland 

35 Fachbereich Physik, Bergische Universitdt Gesamthochschule Wuppertal, Wuppertal, 
Germany"" 

36 DESY, Institut fur Hochenergiephysik, Zeuthen, Germany" 

37 Institut fur Teilchenphysik, ETH, Zurich, Switzerland!" 

38 Physik- Institut der Universitdt Zurich, Zurich, Switzerland} 

39 Visitor from Physicss Dept. University Louisville, USA 

40 On leave from LBL, Berkeley, USA 

41 Institut fur Physik, Humboldt- Universitdt, Berlin, Germany" 

42 Rechenzentrum, Bergische Universitdt Gesamthochschule Wuppertal, Wuppertal, 
Germany" 

43 Physics Department, University of Arizona, Tuscon, USA 
t Deceased 

" Supported by the Bundesministerium fur Bildung, Wissenschaft, Forschung und 
Technology, FRG, under contract numbers 6AC17P, 6AC47P, 6D057I, 6HH17P, 
6HH27I, 6HD17I, 6HD27I, 6KI17P, 6MP17I, and 6WT87P 

b Supported by the UK Particle Physics and Astronomy Research Council, and formerly 
by the UK Science and Engineering Research Council 
c Supported by FNRS-NFWO, IISN-IIKW 

d Supported by the Polish State Committee for Scientific Research, grant nos. 

115/E-743/SPUB/P03/109/95 and 2 P03B 2U 08p01, and Stiftung fur 

Deutsch-Polnische Zusammenarbeit, project no. 506/92 

e Supported in part by USDOE grant DE F603 91ER40674 

? Supported by the Deutsche Forschungsgemeinschaft 

9 Supported by the Swedish Natural Science Research Council 

h Supported by GA CR grant no. 202/93/2423, GA AV CR grant no. 19095 and GA 



3 



UK grant no. 342 

1 Supported by the Swiss National Science Foundation 
i Supported by VEGA SR grant no. 2/1325/96 



4 



1 Introduction 



The electron-proton collider HERA allows to explore a new kinematic region in deep 
inelastic scattering (DIS) down to the very small Bjorken-x of about 10~ 5 , where new 
dynamic features of QCD may show up. So far, the description of the nucleon structure 
function measurement by perturbative QCD, cast into the DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) parton evolution equations |]J, has been extremely successful, 
and constitutes one of the major successes of QCD. At small enough x, however, these 
equations are expected to cease to be a good approximation. For the small x regime the 
BFKL (Balitsky-Fadin-Kuraev-Lipatov) equation |2| has been suggested, and it would be 
very interesting to test QCD in such a new regime. To lowest order the BFKL and DGLAP 
equations resum the leading logarithmic (a s ln(l/x)) n and (a s ln(Q 2 /Qg)) n contributions 
respectively, where Q 2 is the invariant mass squared of the virtual photon and Qq is the 
cut-off for the perturbative evolution. In these approximations the leading diagrams are 
of the ladder type (Fig. p]). The leading log DGLAP resummation corresponds to a strong 
ordering of the transverse momenta k<r (w.r.t. the proton beam) in the parton cascade 
(Qo ^ ^Ti ^ •••&r» ^ ■■■Q 2 )- In the BFKL regime, the transverse momenta follow a kind 
of random walk (k^ ~ k^ i+1 ) |§. 




Figure 1: Generic diagram for parton evolution. 

The strong rise towards small x of the structure function F<i measured at HERA j|, [5[] 
can be described with DGLAP evolution || [7j and is also consistent with the BFKL 
expectation. Less inclusive measurements, namely of the hadronic final state emerging 
from the cascade, may offer more sensitive means of searching for BFKL evolution. One 



5 



possible signature is the quantity of transverse energy, Et, produced in a region between 
the current (struck) quark and the proton remnant ||. As a consequence of the strong 
kr ordering, DGLAP evolution is expected to produce less Et than BFKL evolution §. 
The HERA Et flow data [|1(J can be interpreted consistently within the BFKL framework, 
but from these data alone it is not possible to disentangle non-perturbative hadronization 
effects from perturbative parton radiation, precluding an unambiguous test of the underly- 



ing parton dynamics 12|. In fact, the data can also be described with a DGLAP based 



parton shower model with the assumption that the hadronization effects are large W2 



In this paper complementary measurements of charged particle transverse momentum 
(Pt) spectra are presented. Such spectra have been proposed [13] as a more direct probe 
of the underlying parton dynamics than the Et flow measurements. It has been shown 
with the aid of QCD models that the high-p^ tail is sensitive to parton radiation, whereas 
the production of hard particles from hadronization is suppressed. For the kT unordered 
scenario a harder pt spectrum is expected in the central region than for the fey ordered 
evolution, because parton emissions are less restricted in k T . The unsuppressed gluon 
radiation should give rise to a hard tail of the pt distribution for the hadrons emerging 
from the cascade. The tail would be absent if the Et seen in the data stems predominantly 
from hadronization. A particular good discrimination can be obtained with events in 
which Et is large. Radiation of a hard gluon may lead to large Et and a high pt particle, 
but if hadronization were responsible for the large Et, many soft particles are expected 
instead. Furthermore, a different x-behaviour of the pt spectra is predicted. The phase 
space increases with W 2 ~ Q 2 /x (W is the invariant mass of the hadronic final state) 
and therefore one may expect more gluon radiation, more Et and harder pt spectra 
with decreasing x, but not if kT ordering or another mechanism is invoked to suppress 
gluon radiation. The measurements presented here are compared to predictions from 
Monte Carlo models with different mechanisms for gluon radiation. They incorporate 
QCD evolution in different approximations and utilize phenomeno logical models for the 
non-perturbative hadronization phase. 

The MEPS model (Matrix Element plus Parton Shower) incorporated within the 
LEPTO generator flUl , uses the first order QCD matrix elements, with additional soft 
emissions generated by adding leading log DGLAP parton showers. Though this ap- 
proach is widely used with considerable success, it has to be pointed out that the DGLAP 
formalism was derived for the totally inclusive structure function evolution, and it can 



be questioned whether it is fully applicable for exclusive final states [12|. In the colour 



dipole model (CDM) ]T5|] , as implemented in ARIADNE [IB|, radiation stems from colour 
dipoles formed by the colour charges. Both programs use the Lund string model Wj] for 



hadronizing the partonic final state. The HERWIG model |18[ is also based on leading 
log parton showers, with additional matrix element corrections jn|. This model imple- 
ments an alternative (cluster) fragmentation scheme [pOR . The CDM description of gluon 
emission is similar to that of BFKL evolution to the extent that the gluons emitted by 
the dipoles do not obey strong ordering in fc^ [HI- In the MEPS and HERWIG models 
the partons emitted in the cascade are strongly ordered in fc^, because they are based on 
leading log DGLAP parton showers. The models are used in their versions LEPTO 6.4 
for MEPS, ARIADNE 4.08 for CDM and HERWIG 5.8, together with the parton density 
parameterization GRV . They provide a satisfactory overall description of current DIS 
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final state data [22|, in particular of the Et flows []. For this analysis the use of the 
alternative parameterization MRS-H |24j would give model predictions which differ only 
slightly from those shown. 



2 Detector Description 



A detailed description of the HI apparatus can be found elsewhere [2JJ]. The following 
briefly describes the detector components relevant to this analysis. 

The hadronic energy flow and the scattered positron are measured with a liquid 
argon (LAr) calorimeter and a backward electromagnetic lead-scintillator calorimeter 
(BEMC) respectively. The LAr calorimeter |26| extends over the polar angular range 



4° < 9 < 153° with full azimuthal coverage, where 9 is defined with respect to the pro- 
ton beam direction (+z axis). The depth varies between 4.5 and 8 interaction lengths 
for 9 < 125°. Test beam measurements of the LAr calorimeter modules show an energy 



resolution of &e/E « 0.50/ \J E [GeV] © 0.02 for charged pions p7| . The absolute scale of 
the hadronic energy measurement is presently known to 5%, as determined from studies 
of the transverse momentum (pt) balance in DIS events. 

The BEMC |28[ (with a depth of 22.5 radiation lengths and 1 interaction length) 
covers the backward region, 151° < 9 < 177°. A major task of the BEMC is to trigger 
and to measure precisely scattered positrons in DIS processes with Q 2 values ranging 
from 5 to 100 GeV 2 . The BEMC energy scale for electrons is known to an accuracy 
of 1%. Its resolution is given by a E /E = 0.10/ \Je [GeV] © 0.42/£[GeV] © 0.03 |28[. A 
backward proportional chamber (BPC), in front of the BEMC, with an angular acceptance 
of 155.5° < 9 < 174.5°, serves to identify electrons and to measure precisely their direction. 
Using information from the BPC, the BEMC and the reconstructed event vertex, the polar 
angle of the scattered positron is known to a precision of 1 mrad. 

The calorimeters are surrounded by a superconducting solenoid which provides a uni- 
form magnetic field of 1.15 T parallel to the beam axis in the tracking region. 

The tracking system consists of a central and a forward part. The central tracking 
system is mounted concentrically around the beam line, covering polar angles of 20° < 9 < 
160°. Charges and momenta of charged particles are measured by two coaxial cylindrical 
drift chambers (central jet chamber, CJC), providing up to 56 space points in the radial 
plane. Longitudinal coordinates are obtained via charge division from the CJC, and from 
dedicated drift chambers which are interleaved with the CJC. For the CJC the resolutions 



achieved /Pt ~ 0.009-pr [GeV] ©0.015 and a e = 20 mrad f25[. The central tracking 

system is complemented by the forward tracker, which is built in three sections covering 
a polar angle range of 7° < 9 < 25°. Each section is made up of a series of multiwire 
proportional and drift chambers which are arranged to facilitate track reconstruction at 
small angles to the proton beam. For tracks fitted to a vertex the resolution has been 



1 However, in LEPTO (among other things) the new concept of soft colour interactions Jl^l had to 
be introduced to reproduce the level of Et seen in the data [lj] . Intriguingly, this mechanism also 

produces rapidity gap events Q at a rate of about 10%, comparable to observation Q. Also the cluster 
fragmentation scheme in HERWIG produces rapidity gap events. In this paper rapidity gap events are 
excluded. 
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shown to be cf Pt /pt ~ 0.02 • pr [GeV]/ sm9 © 0.1 and the angular resolution, ag, to be 
better than 1 mrad [pfifl . 

Scintillation counters installed behind the BEMC are used to reject proton induced 
background. Small angle electron/photon taggers are used for luminosity measurements 
and for the study of photoproduction background. 

3 Event and Track Selection 

The data used in this analysis were collected in 1994, with positrons of energy E e = 
27.5 GeV colliding with protons of energy E p = 820 GeV, resulting in a total centre 
of mass energy of ^/s = 300 GeV . The data correspond to an integrated luminosity 
of 1.3 pb _1 . For this analysis DIS events with 5 < Q 2 < 50 GeV 2 are used, in which 
the scattered positron is observed in the BEMC. The events are triggered by requiring a 
cluster of more than 4 GeV in the BEMC. After reconstruction, DIS events are selected 
in the following way: 

• The scattered positron, defined as the most energetic BEMC cluster, must have an 
energy E' e larger than 12 GeV and a polar angle 8 e below 173° in order to ensure 
high trigger efficiency and a small photoproduction background [|J. 

• The lateral size of the positron cluster, calculated as the energy weighted radial 
distance of the cells from the cluster centre, has to be smaller than 4 cm. The 
cluster must be associated with a reconstructed BPC space point Further reduction 
of photoproduction background and the removal of events in which an energetic 
photon is radiated off the incoming positron is achieved by requiring J2j {Ej — p z j) > 
30 GeV [|j], with the sum extending over all particles j (measured calorimetrically) 
in the event. 

• The radial coordinate of the BPC hit must be less than 60 cm, corresponding to a 
positron angle above 157° with respect to the nominal interaction point, ensuring 
full containment of the positron shower in the BEMC. 

• The z position of the event vertex reconstructed from charged tracks has to be 
within 30 cm of the average of all collision events. 

• The energy in the forward region (4.4° < 9 < 15°) has to be larger than 0.5 GeV 
in order to exclude diffractive events with large rapidity gaps in the forward region 

• Remaining background is rejected by requiring no veto from the time-of-ffight coun- 
ters. 

The kinematic variables are determined using information from the scattered positron: 
Q 2 = 4 E e E' e cos 2 (0 e /2) and y = 1 - (E' e /E e ) ■ sin 2 (fl e /2). The scaling variable x is then 
derived via x = Q 2 /{ys), and the hadronic invariant mass squared is W 2 = sy — Q 2 . 



8 



Kin. bin 


x /10~ 3 


Q 2 / GeV 2 


(x) /10~ 3 


(Q 2 )/ GeV 2 


(W 2 )/ GeV 2 


^events 





0.1-10 


5-50 


1.14 


18.3 


24975 


59463 


i 

X 


n 1-0 9 


5-1 n 

O 1U 


n 1 r 


7 n 


459QR 


4853 


9 


n 9-0 5 


fi-i n 


n 9Q 


8 8 


^1fi8fi 


(S9Q4 


3 


0.2-0.5 


10-20 


0.37 


13.1 


36893 


9564 


4 


0.5-0.8 


10-20 


0.64 


14.0 


22401 


7129 


5 


0.8-1.5 


10-20 


1.1 


14.3 


13498 


7964 


6 


1.5-4.0 


10-20 


2.1 


15.3 


7543 


4720 


7 


0.5-1.4 


20-50 


0.93 


28.6 


32390 


8349 


8 


1.4-3.0 


20-50 


2.1 


31.6 


16025 


6752 


9 


3.0-10 


20-50 


4.4 


34.7 


8225 


3838 



Table 1: The kinematic intervals and the average x, Q 2 and W 2 from the raw data. The numbers of 
events surviving the selections are given. 

• As the precision of the y measurement degrades with 1/y, a cut y > 0.05 is imposed. 
Events in which the positron is poorly reconstructed, or in which an energetic photon 
has been radiated from the incoming positron, are removed by demanding that the 
value of y determined using the hadrons y^ = J2(E — p z )/2E e , where the sum runs 
over all hadronic energy deposits, also be greater than 0.05. 

Charged particle tracks are measured in the forward tracker and the CJC. They are 
required to originate from the primary vertex and their polar angle must lie between 8° 
and 155°. The upper limit ensures that the scattered positron does not enter the sample 
of charged hadronic tracks. Central tracks are required to have a radial track length of at 
least 10 cm to ensure good momentum resolution and a measured point in the inner CJC. 
Forward tracks are required to be well reconstructed and to have an acceptable link x 2 in 
the overlap region between the forward and central trackers. After having selected such 
high quality tracks, the efficiency for finding a genuine primary track with p x ^ h > 0.15 
GeV is better than 93% for central tracks and 70% for forward tracks. 

4 Results 

The results are presented for the hadronic centre of mass system (CMS), i.e. the rest 
system of the proton and the exchanged boson, where the direction of the exchanged 
boson defines the positive z' axis. The event sample is divided into 9 different kinematic 
bins, resulting in three slices with almost constant Q 2 and varying x (see Table p. The 
same binning is used for the measurement of transverse energy flows |j0j . 

The data are corrected bin-by-bin for detector effects, including photon conversions, 
secondary interactions, geometrical acceptance and for decay products of particles with a 
lifetime greater than 8 ns (e.g. — > 7r + 7r~ and A — > pir), as well as for QED radiative 
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effects. The corrections are determined using Monte Carlo generated events f\ with a full 
simulation of the HI detector response. The detector simulation has been checked by 
comparing in great detail distributions of track quality properties to the data. A visual 
scan of tracks in real and in simulated events leads to the conclusion that the tracker 
efficiency is simulated to an accuracy of better than 2% for central tracks and 10% for 
forward tracks. No spurious tracks due to noise hits have been found. No significant 
difference was found between the px spectra of positively and negatively charged tracks. 

Transverse momentum spectra are measured with the central tracker for the pseudo- 
rapidity interval 1.5 < rj < 2.5, where rj = — lntan(#*/2), and 9* is the angle with respect 
to the virtual photon direction. The lower rj limit is determined by the requirement of 
good acceptance in all 9 kinematic bins, and the upper limit restricts the measurement to 
a region away from the current fragmentation region. That is where significant differences 
between the different QCD evolution scenarios are expected. 

The measured spectra, steeply falling with px, are shown in Fig. |2| for the nine kine- 
matic bins and for the combined event sample. All distributions shown in this paper are 
normalized to the number of events N that enter the distribution. The inner error bar 
represents the statistical error and the outer error bar the quadratic sum of statistical and 
systematic errors (except for Fig. [|, where only the total error is shown). We first discuss 
the systematic errors and then turn to a discussion of the results. 

The correction factors applied to the original distributions are never larger than 2 or 
smaller than 0.5 and in most cases they are between 0.9 and 1.1. The pt bins are much 
larger than the experimental resolution and bin centre corrections are applied. The boost 
to the CMS introduces a well understood change of the pt of the tracks of less than 10% 
from the laboratory value. The largest corrections are those necessary to compensate 
for event migration effects (contributing mainly at large px) and track selection cuts and 
detector acceptance (contributing mainly at small px)- The main sources of systematic 
uncertainty are: 

• The model dependence of the acceptance corrections (ARIADNE versus LEPTO). 
These differ by as much as 20%, but in general agree to within 5%. 

• The track selection cuts. Variation of these leads to a 10% systematic error for par- 
ticles with px < 0.5 GeV and for the kinematic bins at low W, where measurements 
are close to the limit of the detector acceptance. 

• A contamination of up to 5% from photoproduction for the kinematic bins at large W 
[ |33|| . In order to estimate the uncertainty due to this contamination, stricter electron 
quality cuts have been applied at the cost of statistics, and events with signals in 
the electron/photon taggers are rejected. Differences in the result obtained with 

2 The event generator used is DJANGO |5(J which simulates electroweak interactions using the 
HERACLES Q algorithm and QCD corrections using the ARIADNE [pj program. The detector 
response is simulated using a program based on the GEANT |}2| package. The Monte Carlo events are 
reconstructed in the same manner as the e + p collision data. The simulated and reconstructed events 
obtained using DJANGO describe the raw data spectra very well. Differences in the correction factors 
obtained with the alternative model LEPTO, which describes the data less well, are reflected in the 
systematic errors. 
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respect to the standard selection are reflected in the systematic error and are up to 
30% in kinematic bin 1, below 10% in bin 2 and below 5% in the other bins. 



• The uncertainty of the BEMC energy scale. This affects the boost to the CMS. 
However, the effect is negligible for all bins except for those at large x and Q 2 , 
where it may give rise to an uncertainty up to 15%. 

The pt spectra shown in Fig. ^ are compared to the models ARIADNE, LEPTO 
and HERWIG, representing the suppressed (LEPTO, HERWIG) and unsuppressed (ARI- 
ADNE) gluon radiation scenarios. At "large" x and Q 2 (here "large" means Q 2 ~ 35 GeV 2 
and x ~ 0.004) all models provide a good description of the data. At smaller x and Q 2 , 
LEPTO and HERWIG fall significantly below the data for p T > 1 GeV. ARIADNE gives a 
good description of the data over the full kinematic range. As shown in ||13|| , the hardness 
of the pt spectrum can be related to the parton radiation. The shortfall of the models 
with suppressed gluon radiation indicates that at small x there is more high parton 
radiation present than is produced by the models based upon leading log DGLAP parton 
showers. The data are well described by the ARIADNE model in which parton radiation 
is more abundant. 

In order to look for the predicted hardening of the px spectrum with decreasing x due 
to increasing gluonic activity, the pt spectra are compared in Fig. |3]a for the highest and 
lowest x bins (differing by a factor 6 in x) of the intermediate Q 2 (~ 14 GeV 2 ) slice. 
There is a trend consistent with the prediction of the ARIADNE model, however not very 
significant. Because W 2 ~ Q 2 /x, this effect can also be understood as a.W dependence. 

From a comparison of bins 2 and 7 (Fig. ^|b), where W 2 is fixed (~ 32000GeV 2 ) and 
Q 2 varies by a factor 3 (from 9 to 29 GeV 2 ), it can be seen that the px spectrum becomes 
harder with increasing Q 2 . When x is fixed (x = 0.0021), and Q 2 increased by a factor 
2, the pt spectrum also becomes harder, see Fig. ^|c. We have observed previously that 
in the current fragmentation region the transverse energy, Et, increases with Q 2 for fixed 
W [0 and that the Q 2 dependence diminishes towards the proton remnant direction. The 
Et measured in the interval 1.5 < r] < 2.5 still increases slightly with Q 2 . In the central 
interval —0.5 < i] < 0.5 no Q 2 dependence of is seen for Q 2 < 50 GeV 2 . To summarize, 
at fixed Q 2 there appears to be a hardening of the pr spectrum with decreasing x and 
correspondingly increasing W, while at fixed W there is a hardening for increasing Q 2 . 

Both the forward tracker and the CJC are used to measure the pr spectrum in the 
more central region 0.5 < 77 < 1.5 (see Fig. 4), where the difference between the two types 
of models becomes stronger. Bins 6 and 9 contain roughly equal numbers of forward and 
central tracks whereas the remaining bins are based on central tracks alone since they 
contain few forward tracks. In addition to the sytematic errors described above, the 
data in bins 6 and 9 have a contribution from a 10% uncertainty in the forward tracker 
reconstruction efficiency. The failure of the event generators which are based on leading 
log DGLAP parton showers to describe the hard component present in the data now 
becomes more apparent. ARIADNE describes the data well. Comparing the px spectra 
for the two pseudorapidity intervals in the same kinematic bin (bin 3, low x), it is found 
that the spectrum becomes harder when one moves towards the current direction, see 
Fig. |]d. This is not surprising, as the average k T of the gluons emitted in both scenarios, 
DGLAP and BFKL, increases along the ladder towards the photon vertex Q. Also, 
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radiation associated with the current system, for which the relevant scale is the virtuality 
of the photon (Q 2 ), is expected to contribute more in that direction. 

For the LEPTO predictions in this paper we modified some LEPTO parameters to 
improve somewhat the description of the data over the standard setting. The divergency 
cut-off parameters for the matrix element were changed 0, and a fragmentation parameter 
set tuned to the LEP data [3~o] was used. The most sensitive parameter, however, is 
the virtuality cut-off below which final state radiation is not allowed, and fragmentation 
takes over. If this parameter would be changed from its default value of 1 GeV 2 to 3 
GeV 2 , a much better description of the px spectra for 1.5 < r] < 2.5 could be achieved. 
However, this setting would result in a poor description of the observed scaling violations 



in the Feynman-x spectra || ^1 • m addition, the pr spectra measured more centrally at 
0.5 < T] < 1.5 are still much harder than in the modified model. 

In order to obtain a measure for the hardness of the pt spectrum over the entire range 
of pseudorapidity accessible to this analysis, the average multiplicity of charged particles 
with p T > 1 GeV is measured (see Fig. as a function of pseudorapidity. The surplus 
of hard particles in the data over the DGLAP-like models at small x and away from the 
current region is obvious. The best description of the data is achieved by ARIADNE. In 
contrast to the multiplicity of hard particles, the overall multiplicity, which is dominated 
by soft particles, is much better (although not perfectly) described by all models, see 
Fig. H Only the HERWIG model overshoots the data considerably at small x and towards 
the central region. 

The discrimination between the different scenarios can be enhanced by selecting events 
with substantial hadronic activity because, in such events, there is also substantial par- 
tonic activity for the model with unsuppressed radiation, as opposed to the other models. 
Here the transverse energy measured with the calorimeter in the 77 range from 0.0 to 2.0 
is required to be larger than 6 GeV. The pt distribution of the track with the largest 
transverse momentum in the range 0.5 < r] < 1.5, p™ ax , in each event is then determined, 
see Fig. [7|. In this way one can reveal the correlation between large Et production, which 
could be produced by many soft particles in the hadronization, and high pt particles, 
which indicate hard parton radiation. The events contain much harder particles than are 
produced in the models MEPS and HERWIG. Again, ARIADNE with enhanced parton 
radiation describes the data very well. The discrepancy between the data and ARIADNE 
on the one hand, and LEPTO and HERWIG on the other, becomes larger for smaller 
x and Q 2 . It is noteworthy that the colour dipole model (ARIADNE) is able to describe 
kinematic regions at both large and small x simultaneously. 



5 Conclusions 

Transverse momentum spectra of charged particles in deep inelastic scattering are mea- 
sured as a function of x and Q 2 in the current and the central fragmentation regions. The 
average charged particle multiplicity is also measured as a function of pseudorapidity. The 
data are compared to QCD models incorporating different schemes for parton evolution in 
the proton, namely both with and without transverse momentum ordering of the parton 

3 parl(8) changed from 0.01 to 0.04 and parl(9) from 1 to 4 (3^. 
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emissions. The latter yields a larger amount of high kr parton radiation between the 
current and remnant systems. In this respect the results obtained using the models stud- 
ied can be considered as representative of the expectations for the hadronic final state 
of the DGLAP and BFKL evolution schemes, respectively. At large x, x > 0.001, all 
models provide an overall satisfactory description of the measured pt spectra. At small 
x however and at central rapidity the pt distributions are significantly harder than ex- 
pected from the models with suppressed radiation. The unsuppressed scenario, as realized 
in a colour dipole radiation model, is able to describe the data. The available models for 
the description of the hadronic final state in deep inelastic scattering based on the con- 
ventional leading log DGLAP equations are not able to describe the data. The observed 
hard pr spectra indicate hard parton radiation, which could be provided by BFKL-type 
contributions. Further interpretation of the data has to await firm QCD calculations and 
a Monte Carlo model which explicitly includes the BFKL terms. 
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Figure 2: The transverse momentum spectra of charged particles, measured in the CMS in the 
pseudorapidity interval 1.5 < 77 < 2.5. Data are shown for nine different kinematic bins (see Table 1) 
and the combined sample (bin 0). For comparison, the models ARIADNE (full line), LEPTO (dashed) 
and HERWIG (dotted) are overlayed. The mean values of x and Q 2 are indicated. The inner error bars 
represent the statistical errors, the outer error bars the quadratic sum of statistical and systematic errors. 
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Figure 3: a) Comparison of the pr spectra at high x (bin 6) and at low x (bin 3) for fixed (interme- 
diate) Q 2 . b) Comparison of the pr spectra at two different Q 2 values (bins 2 and 7) for fixed (large) 
W. c) Comparison of the pr spectra at two different Q 2 values (bins 6 and 8) for fixed (large) x. d) 
Comparison of the pt spectra in two different rj intervals in the kinematic bin 3 (low x, intermediate Q 2 ). 
The mean values of x and Q 2 are indicated in the figures. The error bars represent the quadratic sum of 
the statistical and systematic errors. 
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Figure 4: The transverse momentum spectra of charged particles, measured in the CMS in the 
pseudorapidity interval 0.5 < rj < 1.5. Data are shown for seven different kinematic bins (see Table 1) 
plus the combined sample (bin 0). For comparison, the models ARIADNE (full line), LEPTO (dashed) 
and HERWIG (dotted) are overlayed. The mean values of x and Q 2 are indicated. The inner error bars 
represent the statistical errors, the outer error bars the quadratic sum of statistical and systematic errors. 
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Figure 5: The CMS pseudo-rapidity spectra for charged particles with px > 1 GeV. The proton 
remnant direction is to the left. Data are shown for nine different kinematic bins (see Table 1) plus the 
combined kinematic region (bin 0). For comparison, the models ARIADNE (full line), LEPTO (dashed) 
and HERWIG (dotted) are overlayed. The mean values of x and Q 2 are indicated. The inner error bars 
represent the statistical errors, the outer error bars the quadratic sum of statistical and systematic errors. 
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Figure 6: The CMS pseudorapidity spectra for all charged particles. The proton remnant direction is to 
the left. Data are shown for nine different kinematic bins (see Table 1) plus the combined kinematic region 
(bin 0). For comparison, the models ARIADNE (full line), LEPTO (dashed) and HERWIG (dotted) are 
overlayed. The mean values of x and Q 2 are indicated. The inner error bars represent the statistical 
errors, the outer error bars the quadratic sum of statistical and systematic errors. 
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Figure 7: The p™ ax distribution as explained in the text, for the interval 0.5 < rj < 1.5. Data are 
shown for seven different kinematic bins (see Table 1) plus the combined kinematic region (bin 0). For 
comparison, the models ARIADNE (full line), LEPTO (dashed) and HERWIG (dotted) are overlayed. 
The mean values of x and Q 2 are indicated. The inner error bars represent the statistical errors, the 
outer error bars the quadratic sum of statistical and systematic errors. 
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